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Abstract 



o; 

We discuss an embedding of a vector field for the nonholonomic Routh sphere into a 
subgroup of commuting Hamiltonian vector fields on six dimensional phase space. The 
corresponding Poisson brackets are reduced to the canonical Poisson brackets on the Lie 
^3 ■ algebra e*(3). It allows us to relate nonholonomic Routh system with the Hamiltonian 

■ system on cotangent bundle to the sphere with canonical Poisson structure. 

m 

1 Introduction 

Let us consider a smooth manifold M. with coordinates x\,...x m and a dynamical system 
defined by the following equations of motion 

Xi = Xi , i = 1, . . . , m. (1.1) 

We can identify this system of ODE's with vector field 



which is a linear operator on a space of the smooth functions on Ai that encodes the infinitesimal 
evolution of any quantity 

along the solutions of the system of equations (|1.1[) . 

In Hamiltonian mechanics a Hamilton function H on M. generates vector field X describing 
a dynamical system 



X = X H =PdH. (1.3) 

Here dH is a differential of H and P is some bivector on the phase space M. . By adding some 
other assumptions we can prove that P is a Poisson bivector. In fact it is enough to add energy 
conservation 

H = X H (H) = (PdH, dH) = 
and compatibility of dynamical evolutions associated with two functions -ffi,2 

X Hl (X H2 (F)) = X H2 (X Hl (F)) + X Xhi (H2) (F) , 

see [U [13] and references within. 

For a considerable collection of nonholonomic dynamical systems vector fields are created 
using Hamilton function H and density g of the invariant with respect to X measure 

X = gPdH . (1.4) 
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This vector field X (|1.4j) is a so-called conformally Hamiltonian vector field, see examples of 
such fields in HH H [501 [21]. 

Below we discuss nonholonomic Routh sphere problem with the vector field X, which is a 
sum of commuting vector fields PdHk determined by integrals of motion H±, . . . , H n : 

X = g 1 PdH l + ■■■+ g n PdH n . (f .5) 

Such decompositions are well-known in bi-Hamiltonian geometry |18j and for some other non- 
holonomic dynamical systems [5J [2U] . If we have decomposition (II .5[) we can say that common 
levels of integrals of motion form a Lagrangian foliation associated with the Poisson bivector 
P with all the ensuing consequences. We want to highlight that we discuss only properties of 
foliations and do not discuss linearization of the corresponding flows. 

It is known that a proper momentum mapping for the non-Hamiltonian vector field X 
associated with the Routh sphere has a " focus- focus" singularity [TU]. According to [TT] the 
nontriviality of the corresponding monodromy is a coarest obstruction to the existance of global 
action-angle variables. 

2 The Routh sphere 

Following to [H [HI EH E3 let us consider a rolling of dynamically symmetric and non balanced 
spherical rigid body, the so-called Routh sphere, over a horizontal plane without slipping. 
Dynamically non balanced means that that the geometric center differs on the center of mass, 
whereas dynamically symmetric means that two momenta of inertia coincide to each other, for 
instance I\ = I 2 - The line joining the center of mass and the geometric center is axe of inertial 
symmetry. 

The moving sphere is subject to two kinds of constraint: a holonomic constraint of moving 
over of a horisontal plane and no slip nonholonomic constraint associated with the zero velocity 
in the point of contact 

v + ioxr = 0. (2.1) 

Here uj and v are the angular velocity and velocity of the center of mass of the ball, r is the 
vector joining the center of mass with the contact point and x means the vector product in IR 3 . 
All the vectors are expressed in the so-called body frame, which is firmly attached to the ball, 
its origin is located at the center of mass of the body, and its axes coincide with the principal 
inertia axes of the body. 

In the body frame the angular momentum M of the ball with respect to the contact point 
is equal to 

M = Iqw , I Q = I + mr 2 E -mr®r. (2.2) 

Here E is a unit matrix, m is a mass and I = diag(ii, 12,1s) is an inertia tensor of the rolling 
ball, respectively. 

If 7 = (71, 72, 73) is the unit normal vector to the plane at the contact point, then 

r = (Ryi, Rj2, -R73 + a) , 

where R is a radius of the ball and a is a distance from geometric center to the center of mass. 

The phase space, initial equations of motion, reduction of symmetries and deriving of the 
reduced equation of motion are discussed in [H [3 [TUJ. We omit this step and begin 
directly with the reduced equations of motion on the six dimensional phase space M. with local 
coordinates x = (7, M) : 

M = M x lj + mr x(wxr), 7 = 7x0;. (2.3) 
A straightforward calculation shows that these equations (|2.3p possess four integrals of motion 
Hi = (M, uj) , H 2 = (M,M)-mr 2 H 1 , H 3 = (M,r) ^ = (7,7), (2.4) 
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and the following invariant measure 

It = g-\i) d 7 dM, 5(7) = yjhh + hmB?{Hi - 7 |) + hm(R l3 + a) 2 • (2.5) 

As for the symmetry Lagrange top there are two linear in momenta integrals of motion. First 
integral 

H 3 = {M,r) 

is a well-known Jellet integral [12| . see also §243, p. 192 in the Routh book [T5] . Second 
integral 

5(7)143, (2.6) 

was found by Routh in 1884 [15] and recovered later by Chaplygin in 9 . These linear integrals 
are related with quadratic in momenta integral of motion 

Ho = Ji H\ 4- — Ho — — Ho . 

h 2 h 

According to 4] namely this quadratic integral survives at I\ ^ I 2 in contrast with the linear 
integral Hi. 

If the Routh sphere rolls on a horizontal plane under the influence of a constant vertical 
gravitational force, equations of motion read as 

dU 

M = M x ui + mr x (oj x r) + 7 x — — , 7 = 7x11). (2.7) 

07 

Here U = — mg(r, 7) and g is a gravitational acceleration. These equations possess the same 
invariant measure [i and integrals H 3 , H 4 ,H 2 , whereas quadratic integrals of motion are equal 
to 

tfi = [w,M)+2U, Hi = M 2 -mr 2 (uj,M) -2mg(aIi l3 -mR 2 j 2 (a 2 +j 2 R 2 )y (2.8) 

In this case vector field X defined by equations (|2.7[) has homoclinic trajectories when the Routh 
sphere is either spinning very slowly about a vertical axes which passes through the center of 
mass and the geometric center, either axis is in the same position but the value of the Jellet 
integral is slightly less than the threshold value [IS] . 

2.1 Poisson brackets 

For the Routh sphere six equations of motion (|2.3I) possess four integrals of motion and an 
invariant measure and, therefore, they are integrable by quadratures according to the Euler- 
Jacobi theorem. It allows us to suppose that common level surfaces of integrals form a direct 
sum of symplectic and lagrangian foliations of dual dynamical system which is hamiltonian with 
respect to the Poisson bivector P, so that 

[P, P] = , PdC 1>2 = PdHij = , (PdH h dH k ) = {Ht, H k } = . (2.9) 

Here [.,.] is the Schouten bracket and (i,j,l,m) is the arbitrary permutation of (1,2,3,4). In 
fact, here we suppose that the Euler-Jacobi integrability of non-Hamiltonian system (I2.3[) is 
equivalent to the Liouville integrability of the dual Hamiltonian dynamical system with the 
same integrals of motion, see [18] . 

The first equation in ()2.9[) guaranties that P is a Poisson bivector. In the second equation 
we define two Casimir elements Hi and Hj of P and assume that rankP = 4. It is a necessary 
condition because by fixing its values one gets the four dimensional symplectic phase space of 
our dynamical system. The third equation provides that the two remaining integrals Hi and 
Hk are in involution with respect to the Poisson bracket associated with P. 

In this note we discuss solutions of the equations (I2.9|) in the space of the linear in momenta 
Mi bivectors P at the different choice of the Casimir functions: 
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1. Ci = (7,7), C 2 = (M,r), He = Hi , H k = H 2 ; 



2. d = (7, 7) , C 2 = g{l) uj 3 Hi = H\, H k = H 3 ; 

3. d = (M, r) , C 2 = 3(7) w 3 H t = H lt H k = H 4 . 
In generic case linear in momenta Casimir functions look like 

Ci,2 = ai,2 (7)7) + &i,a(M,r) +01,25(7)^3 , 121,2,61,2^1,2 e C. (2-10) 

However, the corresponding complete solutions of (|2 .9[) have the same properties as particular 
solutions obtained in the listed above three special cases. 

In [T7] we have solved the same system of equations (|2.9p for the symmetric Lagrange top 
and proved that solutions may be useful for construction of the variables of separation and the 
recursion Lenard-Magri relations for this Hamiltonian system. 

If we have some solution P of (|2.9p . we can get decomposition of initial vector field X by 
commuting Hamiltonian vector fields PdHi and PdH k . The existence of such decomposition 
by the basis of the Hamiltonian vector fields requires to impose one more condition rankP=4. 



3 First Poisson bracket 



Substituting linear in momenta Mj anzats for entries of the Poisson bivector 
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fc=i 



into (|23|) at 



Ci = if 4 = (7,7) , C 2 = H 3 = (M,r) 



Hi = Hi, H k = H 2 



and solving the resulting system of algcbro-diffcrcntial equations one gets the following Propo- 
sition. 

Proposition 1 In this case generic solution of \2.9\) is parameterized by two functions 0(71/72) 
and /3(7 3 ) 

I> 



P = ag[ -rj M a 



(3.2) 



Here matrices ^ a ,/3 o^e equal to 

7i72(fl73+a) 7lCg73+a) 



H(7?+7l) 



R(7?+7l) 



"72 



7i Cg73+a) 7l72(-R73+a) 





7172 73 7i73 n \ 

72 73 7172 73 Q 



72 



-71 



0/ 



and sfcew symmetric matrices M a .p have the form 
M a = 



/ Q (71^+72^2)^73 + 0) 



^(7i 2 + 7 2 2 ) 




-Af a 
Mi 
/ 



M, 



M 3 



73(71 Mi + 72M2) mo-R(R-f 3 + a) ma^ 2 R 2 \ 



7?+7 2 2 



V 



3 Z 




where g = g("f) and 



a = mR 



(m(r, j)C 2 + 13(71 Mi + 72M2) + /1M373J 
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The proof is a straightforward solution of (|2.9j) using linear in momenta anzats. 
The corresponding Poisson brackets read as 

«ff7i 72(^73 +a) , /?7i7273 r ,, t «57i(^73 + a) , 



|iWl,7i} = — g— j— 1 2") {-lWl,72}— / 2 , 2\d 1 2 "■ 

(7i+7 2 )^ 7i+7 2 (7i+72)^ 7i+7 2 

{Mi , 73 } = -/37 2 , { M 2 , 73 } = /?7i , 

r,r 1 _ Ugij{Rj3 + a) Pill* r M l _ «g7l 72(-R73 + O) /?7l7273 

1 2,711 ~ (7! 2 + 7 2 2 )i? 7 1 2 + 7l' ^,72}- (7 2 +7 2 )i? 7 2 + 7r 

{M 3 ,7i} = «57 2 , {M 3 ,7 2 } = -«37i , {M 3 ,7 3 } = 0, {7i,7j} = °) 



(3.3) 



and 



fir 1 a.g(7i^i + 72M 2 )(i?7 3 + a) / 73 (71 Mi + 72 M 2 ) a (#73 + a) 

t™b* J 2} — 7~2~, 2v5 1 3 J - ] 2 2 

(7l+7 2 )^ V 71+72 9 

{Mi, M 3 } = -agM 2 + £-=- 72 , {M 2) M 3 } = a<?Mi - 71 . 

3 5 

In generic case rankP = 4, however if a = or /3 = there are additional Casimir function 73 

and 71/72, respectively. A particular form of this brackets was obtained in [14] , 

Using this Poisson bivector we can obtain basis of the commuting Hamiltonian vector fields 

Xi = PdH x and X 2 = PdH 2 , 

and try to expand initial non-Hamiltonian vector field X (12. 7p by these vector fields. 

Proposition 2 Using Poisson brackets 13. 3\) we can rewrite reduced equations of motion for 
the Routh sphere \2. 7| ) in the following form 

x k = gx{x k ,H\} + g 2 {x k ,H 2 }, k = l,..6, (3.4) 

if and only if 

"(71/72) = const , 0(73) = ag ^1 + -^-^ . (3.5) 

In this case coefficients are equal to 

(Rj 3 + a)/i - -R73/3 a 

91 ~ 2ag(h - / 3 )(i? 7 3 + a) ' 92 ~ 2ag(h - I 3 )(R l3 + a) ' 

The proof is a straightforward verification of the equations (|3.4[) . 

There are other special values of the functions a and f3 according to the following 

Proposition 3 If 

( 1 \ +. a at \ ag(hRj 3 - I 3 {R-f 3 + a) 

0(71/72) = const, and /3( 73 ) = — 7— — ■ rr , (3.6) 

R\{h ~ hm + am(r, 7) ) 

the Poisson bivector P V3.2\l is compatible with the canonical Poisson bivector Pq on the Lie 
algebra e*(3) 

r 
r T m 



Pi 

where 






73 


-72 


-73 





7i 


72 


-7i 






M 3 -M 2 
M I -M 3 Mi 
-Mi 
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The proof is a calculation of the Schouten bracket [Pq, P] =0. 

Remind, compatibility means that a linear combination of this bivectors 

P x = P Q + XP , A 6 C , 

is a Poisson bivector at any value of A. It also means that P (|3.2p is a trivial deformationof Pq, 
see details in [15] . 

It is interesting that namely this condition of compatibility allows us to expand initial 
vector field (|2.7[) by a basis of Hamiltonian vector fields 

±k = gx{xk,Hx} + g 2 {xk,H 2 }, fc=l,..6, (3.7) 

associated with the linear in momenta Routh integral H 2 (|2.6I) . Here coefficients 



1 „ a ( 73(71^1+72^2) 
91 = > 32 = ; 77 ~ TTT, ; — sT I M 3 9 T ^2 



2/3 ' a/j (7^73 - J 3 (#73 + 0)) V 7i + 7 2 2 

depend on coordinates and momenta in contrast with the previous decomposition. 

Summing up, for the Routh sphere integrals of motion are in involution at any values of a 
and f5 

{H 1 ,H 2 } = {H 1 ,H 2 } = 0, Va,j8. 

However, initial vector field X is decomposed by the corresponding commuting Hamiltonian 
vector fields 

X = gi PdHi + g 2 PdH 2 , or X = foPdHx + g 2 PdH 2 



only at the special values of these functions (|3.5j) or (|3.6j) . respectively. 
3.1 Properties of the first Poisson brackets 

Similar to Chaplygin sphere problem |20j and to nonholonomic Veselova problem [21) we can 
reduce this Poisson bracket to the canonical Poisson brackets on the Lie algebra e*(3) and 
identify the Routh sphere model with the Hamiltonian system on two-dimensional sphere. 
One of the possible reductions is given by the following Proposition. 

Proposition 4 After a change of momenta 

T _ 1 ( 1il3(R(li^h+l-2M 2 )-bI- 1 (I l +m(R l3 +a) 2 ) 72 ( 72 Mi- 7 iM 2 ) . 
Ll ~ ^P^l \ ag (R l3 +a) + p + C Tl 

7273(fl(7iMi+72M 2 )-b7 1 - 1 (J 1 +m(fl 7 3+a) 2 ) _ 7l ( 72 Mi- 7l M 2 ) \ (3.8) 



ag(R*13+a) 

^3 = ^+ bm( "l +a) , b = (M,r), c=(L, 7 ), 



C72 



ag agh 

the Poisson brackets {., .} i3. 3\) coincide with the canonical Poisson brackets on the Lie algebra 

e*(3) 

{Li , Lj } Q = e ljk L k , {Li , 7j } Q = e ljk -/ k , {7,; ,7, } Q = , (3.9) 
where Eij k is a completely antisymmetric tensor. 

It is easy to see that the Poisson map (|3.8p is locally defined in the region 

7l 2 + 72 = 1 - ll + 0. 

Namely in this region of the phase space the vector field for the Routh sphere X (|2.7p does not 
have homoclinic orbits HU1. 
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If c = (7, L) = and condition (|3.5p holds, then images of the initial integrals of motion 
are nonhomogeneous second order polynomials in momenta 

a 2 ( g 2 ~hh 2 ff 2 (^73 + a) 2 (£i72 - L 2 ji) 2 
1 7x 2 + 7 2 2 V i? 2 7 3 2 (/i + mr 2 ) 

2bag(R-/ 3 + g)L 3 fr 2 (/i + m{R-f 3 + a) 2 ) 
/ 1 i? 2 ( 7l 2 + 7 2 2 ) + I 2 i? 2 (7 2 + 7 2 ) ' 



(3.10) 



Tr ,' hhr 2 t2 , 5 2 A(i?73 + a) 2 (^i72 -^27i) 2 

til — — n n ^ Li- 



7? + 7l V # 2 ^ 2 7l(^i + mr 2 ) 

2bag(Rj 3 + a)L 3 ( Ii + mr 2 i 2m 



i? 2 (7i 2 +7 2 2 ) ' V^i? 2 (7i 2 + 7 2 2 ) A 

These integrals define a Hamiltonian system on cotangent bundle T*S 2 to the sphere S 2 . 
As for Lagrange top the existence of the linear in momenta integral of motion (|2.6[) 

H 2 = aI\L 3 , 

allows us to explicitly integrate the corresponding Hamiltonian equations of motion by quadra- 
tures. Namely, let us introduce spherical coordinates on the sphere 

sin^cos^ 

7 1 =sm0sin.0, Li = — — p^ - cos <p p g , 

... r cos0cos6» . (3.11) 

72 = cos0sm0, L 2 = — ; — pd, + sm(pp0 , v ' 

sin 9 

73 = cos 6 , L 3 = , 

where <j>, 9 are the Euler angles, p^ and pe are the canonically conjugated momenta 

{4>, P< f,} = {9 lPe } = 1 , {0, 0} = {0, Pe } = {6,p^} = . 
In this variables initial integrals of motion are equal to 

Hi = A{8)p 2 + B{9)p 2 + bC(9) P<p + b 2 V(9), H 2 = -aI lP ^, (3.12) 
where b is a value of the Jellet integral, A, B, C and V are function on 8: 

AOS) =a 2 (l 1+ Ma 3 + 2aficoB fl+ flW 9 ) \ = ff 2 

V / V ^ s / ii+m(a J +2a/t cos #-t-K 2 ) ' 

2ag(R cos l/f^ Ji +m(a 2 +2a.R cos 0+i? 2 ) 

>■ ' /1 fl 2 sin 2 6 ' > i / 2 7?, 2 sin 2 ' 

Thus, similar to the Lagrange top, solution of the Hamiltonian equations of motion generated 
by the Hamilton function and canonical Poisson brackets is reduced to solution of one equation 

9 = 2B(8) Pe = 2^B{8){E - A(8)c 2 - bcC{8) - b 2 V{9)) , 

where E — Hi and c = —H 2 jali are constants of motion. 

Summing up, after reduction of the Poisson brackets to canonical ones we can easily find 
trajectories of the commuting Hamiltonian vector fields 

Xi = PdHi, X 2 = PdH 2 and X 2 = PdH 2 . 

However, in order to get solutions of the Routh equations of motion we have to integrate their 
linear combinations 

X = giXi + g 2 X 2 or X = g t Xi + g 2 X 2 , 

associated with conditions (|3.4I) or (|3.4p . respectively. Solution of this problem remains open 
as of yet. 
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3.2 Conformally Hamiltonian equations of motion 

If the Jellet integral is equal to zero C2 = (M, r) = 0, i.e. if 6 = 0, then integrals of motion Hi 2 
(|3.10[) become homogeneous quadratic polynomials in momenta. In this case we can easily find 
variables of separation gi 2 in the corresponding Hamilton- Jacobi equation, if we diagonalize 
simultaneously two quadratic forms Hi 2 (I3.10[) . Then, using these variables of separation, we 
can rewrite initial vector field X (ll.5[) in the conformally Hamiltonian form (ll.4[) . 

Namely, at C 2 = b = integrals of motion Hi <2 satisfy to the following separated relations 

$ i (q i ,p i ,H 1 ,H 2 )=0, k = l,2. 

Here qi^ and pi.2 are canonically conjugated variables of separation. In this case, according to 
[16] , these integrals Hip are in involution 

{H 1 ,H 2 } f = 

with respect to the Poisson brackets 

Ui,Pi}f = fi(qi,Pi) , {<72,P2}/ = /a (32,^2) , {?i 5 92}/ = {Pi,P2}f = 0, 

labelled by two arbitrary functions fi t2 . The corresponding Poisson bivector Pt and integrals 
of motion -Hi,2 satisfy to the equations 

P f dH t = F a PdHi + F i2 PdH 2 , i = 1,2. (3.13) 

Here functions Fij depend on fx >2 and form the so-called control matrix [T6j [18] . 

Proposition 5 If X hl.5\) is a linear combination of the commuting Hamiltonian vector fields 

X = giPdHi + g 2 PdH 2 , 

and coefficients gi^ are special combinations of Fij 

gt = gfaiFu + a 2 -F i2 ) , i = 1,2, 

there is a Poisson bivector Pf, which allows us to rewrite X in the conformally Hamiltonian 
form 

X = giPdHi + g 2 PdH 2 = gP/dH , H = a x Hi + a 2 H 2 . 
In this case H is a sum of initial physical integrals of motion Hi 2- 

For the Routh sphere at C 2 = variables of separation q\ >2 are functions only on coordinates 
7j. Thus, the desired bivector Pf may be directly obtained from P (13.21) at 

a - -H , p — - - — -— , 
(7,Ir) 

so we have 

X = giPdHi + g 2 PdH 2 = ~ PfdHi . 

zp 

At C2 = variables of separation qi^ have to be functions on coordinates 7$ and momenta Mj 
and, therefore, entries of Pf have to be more complicated functions on Mi. Unfortunately, we 
do not know how to get variables of separation for the nonhomogeneous polynomial integrals 
of motion (j3. 10[) on the sphere. 
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4 Second and third Poisson brackets 

Let us substitute linear in momenta M{ anzats (|3.1|) into the equations (|2.9[) at 

C 1 = J ff 4 = (7,7), 0> = #2 =s(7)w 3> H e = Hi, H k = H 3 



Proposition 6 In this case integrals of motion are in involution {Hg,Hk} = if and only if 
bivector 

zs labelled by three arbitrary functions 0^(7) entering into the matrix 

-72(«1 + 73« 3 ) "72^2 + 7371^3 jl+m^+aV 

mfl7i(7iai+7 2 a2)(fl73+n) 

7l«l 7l"2 /l+m(jR73+a) 2 

7l72«3 ~7l"3 

and into i/ie sfcew symmetric matrix 

M' 12 = -aiA/i - a 2 Af 2 + a 3 (7iM 3 - 73 M x ) 

^371^73 + a)(m 2 (7,r)iJ 3 + m(J 3 (7l M l + 72 Af 2 ) + /173M3 



,9 2 



, mR(Rj 3 + a)(7iai + 7 2 a 2 ) 

h + m(i?7 3 + ay 



a 3 7i72i? 2 (m 2 ( 7 , r)i/ 3 +m(/ 3 (7i Mi + 72 Af 2 ) + /173M3)) 



<J 2 



, uax(ix-j 3 -t- un/yiui -t- 72 u; 2j ; 

IVlo o — — — Ml 

2 ' s t, -x- m rp-v„ _i_ ^2 1 



mR(Rj 3 + a)(7i«i + 72^2) 
h + m(i?7 3 + a) 

a 3 7 2 i? 2 (m 2 ( 7 ,r)iJ 3 + m(l 3 (^ 1 M 1 + 72M 2 ) + /173M3) 

? ' 

If we impose additional restriction rankP' = 4, the third equation [P' , P'\ — in h2. 9\) has the 
following single solution 

h + m(i? 73 + a) 2 llhR(h + ma{R l3 + a)) 

"1 1 7f \ , 

7i 7H I 7)H 

72/1 -R(/i + ma(Rj 3 + a)) 

Q 2 = = \ 1 (4.2) 

(l7,r) 

/l + m{R l3 + a) 2 (7? + 7 2 2 )^i (h + ma{R l3 + a)) 

7173 7173(17. r ) 

The proof is a straightforward solution of the differential equations using substitution (|3.1|) . 

Using this Poisson bivector we can get a basis of Hamiltonian vector fields and an expansion 
of the initial vector field 

X = g[P'dH 1 + g' 2 P'dHl (4.3) 
by these vector fields. The corresponding coefficients are equal to 

, _ f , _ I\{I\ + ma(R-f 3 + a)) - mR(-f, Ir) 

- 9l = ~2 7l a 3 V 92 ~ 2 5 2 i?/ 2 (i?73 + a) ' 
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Similar to the first Poisson bracket (|3.3|) . there is a transformation of momenta Mj — > Li, which 
reduces this Poisson brackets {.,.}' to canonical Poisson brackets on the Lie algebra e*(3) . 
Now let us consider the third possible choice of the linear in momenta Casimir functions 

d = (M,r), C 2 =g(~/)u 3 H e = Ht, H k = H 4 . 

In this case generic solution of (|2.9|) coincides with the previous solution P' (14. ip 

P" = P'\ 

I7iai+7 2 a2=0 ' 

at jioti + 72 «2 = 0. It is easy to see that rankP" = 3 and there is a third Casimir function 

Hi = {l, l) 

P"dC 3 = , C 3 = H 4 . 

Consequently, in this case we have only one nontrivial Hamiltonian vector field P"dHi, which 
does not form a basis. 

In the similar manner we can consider generic linear in momenta Casimir functions (|2.10[) . 

Proposition 7 We can not rewrite equations of motion \2. 7p on the six- dimensional phase 
space in the conformally Hamiltonian form 

X = gPdF(H 1 ,H 2 ,H 3 ,H i ) (4.4) 

using linear in momenta Poisson bivector P satisfying equations H2.9]) . Here F(Hi, H 2 , H 3 , H 4 ) 
is an arbitrary function on integrals of motion for the Routh sphere. 

The proof is a straightforward verification of the fact that common system of equations (|4.4I 
and (|2.9p is inconsistent at the space of the no more then linear in momenta Casimir functions 
and Poisson bivectors. 

5 Conclusion 

It is well known that equations of motion for the nonholonomic Routh sphere are integrable 
by quadratures according to the Euler-Jacobi theorem. We identify the corresponding level 
sets of integrals of motion with the Lagrangian foliations associated with two different Poisson 
bivectors P and P'. The corresponding expansions of the initial vector field X (|3.4[) and (|4. 3[) 

X = giPdH x + g 2 PdH 2 = g[P ' dH x + g' 2 P ' dH\ 

may be considered as a counterpart of the standard Lenard-Magri recurrence relations 

X = PdH x = fiP'dH 2 + f 2 P'dH 3 

for two dimensional bi-Hamiltonian systems (/1 = 1, f 2 — 0), quasi bi-Hamiltonian systems 
(/2 = 0) or bi-integrable systems (V/i^), which appear in Hamiltonian mechanics [TBI [T7l H8] . 

We would like to thank A.V. Bolsinov, A.V. Borisov and I.S. Mamaev for useful discussion 
of applications of the Poisson geometry to the different nonholonomic systems. 
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